Low-temperature magnetization in geometrically frustrated Tb 2 Ti 2 07 
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The nature of the low temperature ground state of the pyrochlore compound Tb2Ti2C>7 remains 
a puzzling issue. Dynamic fluctuations and short-range correlations persist down to 50 mK, as ev- 
idenced by microscopic probes. In parallel, magnetization measurements show irreversibilities and 
glassy behavior below 200 mK. We have performed magnetization and AC susceptibility measure- 
ments on four single crystals down to 57 mK. We did not observe a clear plateau in the magnetization 
as a function of field along the [111] direction, as suggested by the quantum spin ice model. In ad- 
dition to a freezing around 200 mK, slow dynamics are observed in the AC susceptibility up to 4 K. 
The overall frequency dependence cannot be described by a canonical spin-glass behavior. 



PACS numbers: 75.40.Gb; 75.60.Ej; 75.30.Kz 

Geometrical frustration studies in magnetism have 
revealed a large variety of new and exotic magnetic 
phases^. Among the frustrated materials, the rare-earth 
pyrochlore oxides, in which the magnetic atoms are situ- 
ated on a corner-sharing tetrahedra lattice, form a fam- 
ily in which different kinds of anisotropy and / or mag- 
netic exchange can be obtained, thus providing a play- 
ground for studying the role of these parameters in the 
frustration^. Their chemical formula is R2M2O7 where 
R is a rare-earth and M a non magnetic metallic element. 
Among these compounds, Tb 2 Ti 2 07 has been the focus 
of much interes t si nce no magnetic ordering is observed 
down to 50 mKP^I in spite of antiferromagnetic interac- 
tions leading to a Curie- Weiss temperature of -13 K^H In 
this compound, the magnetic Tb 3+ ions have a strong 
Ising-like anisotropy along the local [111] axis. Instead 
of a long range order predicted to occur around 1 KP3, 
short range correlations are observed down to the lowest 
temperature^!. 

It has been proposed that the ground state and first ex- 
cited state of the Tb 3+ ion are two doublets, separated by 
a gap A ~ 18 K 5 9 . In this framework, it was suggested 
that quantum induced crystal field excitations could help 
explain the absence of long-range order in Tt^T^O-^l. 
Thus this compound was considered to be a "quantum 
spin ice" in analogy to dipolar spin ices which are charac- 
terized by a degenerate ground state in which two spins 
point in and two spins point out of a tetrahedron^. It 
follows that the magnetization curves at low enough tem- 
peratures should show a magnetization plateau in the 
presence of an applied field along the [111] direction^ 
as observed in dipolar spin ice compounds^. Notwith- 
standing, recently it has been proposed that a tetragonal 
distortion which occurs at low temperature^ED lifts the 
degeneracy of t he low est doublet, thus generating a sin- 
glet ground stat d 17 * 1 ^ . In this case, the resulting Tb mag- 
netic moment is induced by exchange interactions^. As 
the nature of the Tb 3+ ion ground state is still a debated 
question^, magnetization curves at very low tempera- 
ture could be a key to resolving it. 

Another puzzling feature of the Tb 2 Ti207 low tem- 
perature magnetic state is the apparent contradiction 



between the persistent dynamics down to 50 mKP, and 
the evidence of f reezing observed near 200 mK by AC 
susceptibility^^ and in Zero Field Cooled-Field Cooled 
(ZFC-FC) magnetization measurement*^! which was at 
first attributed to the presence of defects^. Irreversibil- 
ities in muon spin relaxation (/zSR) measurements and 
a minimum in specific heat have been reported at the 
same temperature^. Furthermore, a saturation of the 
muon relaxation rate has been observed below about 2 
KPUD These observations suggest the coexistence of dif- 
ferent time and temperature scales for the spin dynam- 
ics, and up to now, no clear unified picture has emerged 
to describe the dynamics of Tb2Ti207. Although mi- 
croscopic probes (neutron scattering and /zSR measure- 
ments) have been extensively used, detailed studies of the 
magnetization^ and of the AC susceptibility need to be 
performed at very low temperature. 

In this communication we present a systematic study 
of the macroscopic magnetic properties down to 57 mK 
made on Tb2Ti207 single crystals. We show that down to 
this temperature, there is no evidence of a magnetization 
plateau in the [111] direct ion^. We also show that the 
reported glassy phase^H is not a canonical spin glass 
phase, and is strongly dependent on the applied field. 

Magnetization and AC susceptibility measurements 
were performed by the extraction method, using two su- 
perconducting quantum interference device magnetome- 
ters equipped with miniature dilution refrigerators and 
developed at the Institut Neel 25 . AC measurements were 
carried out for frequencies between 1.1 mHz and 2.11 
kHz, with an applied AC field of 0.5 and 4 Oe (No de- 
pendence on the amplitude of the AC field was observed) . 

Three crystals, denoted as B, C, and D, were prepared 
for these measurements (see Suppl. Mat.). These crystals 
have been characterized by specific heat measurements. 
The data for crystal C have been published in Ref. [23] 
The specific heat behaviors of samples B and C are qual- 
itatively the same, with no anomalies, although sample 
B has a slightly larger residual entropy. In contrast, the 
specific heat of sample D is very different with a narrow 
peak around 400 mK, reminiscent of the data published 
in Ref. [2T] These crystals were cut into flat, disk shaped 
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FIG. 1. (Color online) First magnetization curves M vs. fioH 
along the [111] direction between 57 mK and 4.2 K for sample 
B (a) up to 0.5 T, (b) up to 8 T. (c) Comparison between the 
experimental curves at 57 and 100 mK to those predicted in 
Ref. M 



samples of about 30 mg. Four samples were measured: 
two from crystal C, called CI and C2 measured along the 
[111] and [110] directions respectively, one from crystal B, 
and one from crystal D, both measured along the [111] 
direction. Demagnetizing effects are very small due to 
the geometry of the samples (the field orientation was in 
the plane of each disk, except for sample D). The crystal 
alignment with respect to the field is accurate to within 
a few degrees. 

Measurements of the first magnetization curves along 
the [111] direction were carried out on samples B, CI and 
D. Figure [ljb) shows these isotherms M vs. /j-oH for 
fields up to 8 T between 57 mK and 4.2 K for sample B. 
The magnetization reaches 5 ^ts/Tb at 8 T, which is con- 
sistent with previous measurements^!, and corresponds to 
the predicted ground state magnetic moment for the Tb 
iorP. The results are very similar for the three samples 
(See Suppl. Mat.). 

Figure [l|a) shows the magnetization below 0.5 T in 
the same temperature range. These curves can be com- 
pared with the predictions of Molavian and Gingraa^ 
who proposed the emergence of a magnetization plateau 
below 100 mK in the [111] direction. At 57 mK, no clear 
evidence of such a plateau can be observed, which is con- 
firmed by plotting the calculated derivative dM/dH vs 
H (See Suppl. Mat.). Nevertheless a small break of slope 
occurs around /loHt = 0.13 T. This field can be compared 



to the level crossing field £ c =0.082 T in Ref. HI which 
also gives rise to a change of curvature in the calculated 
100 mK magnetization curve (See Fig. [ljc)). In addi- 
tion, we have measured the magnetization curves along 
the [110] direction down to 70 mK in sample C2: they 
are only slightly different from the magnetization curves 
measured along [111] (see Suppl. Mat.). 

Two main differences emerge when we compare our 
measured curves along [111] to the predicted magneti- 
zation curves of Ref. calculated in the independent 
tetrahedra approximation (ITA) (See figurejljc)): i) The 
measured value of the magnetization at 0.5 T is signifi- 
cantly lower than the calculated magnetization, ii) More 
important, our 57 and 100 mK curves are almost iden- 
tical whereas the calculated 50 and 100 mK curves have 
very different shapes. 

Several hypotheses can be proposed to account for 
these discrepancies, i) It has be en sug gested that the ex- 
change interaction is anisotropic! 19 * 26 * In this scheme the 
proposed exchange energies are larger than the isotropic 
exchange constant used in the calculation of Ref. [TU 
Ref. [T2] shows that the increase of the exchange tends 
to decrease the magnetization value at low field and to 
decrease the temperature at which the plateau is formed. 
Thus in this case, our lowest measuring temperature 
could be higher than the temperature for the formation 
of the plateau, ii) It has been recently proposed that 
the Tb ground state is a singletP^^ rather than a dou- 
blet. This scenario presumably precludes the existence 
of the plateau. The computation of the magnetization 
curves within this crystal field scheme would be of great 
interest to compare to the present data, iii) In Ref. [T2l 
the calculation were made in the independent tetrahedra 
approximation. It would be interesting to check the ro- 
bustness of the plateau when including further neighbors. 

We measured the ZFC-FC magnetization as a func- 
tion of temperature for the four samples. All of the data 
show an irreversibility below 250 mK (See inset of Fig. 
[2J: the ZFC curve has a maximum between 150 and 200 
mK, whereas the FC magnetization tends to saturate. 
This behavior is in agreement with previous magnetiza- 
tion measurements^!. It occurs at almost the same tem- 
perature as the irreversibility observed in transverse /iSR 
measurements (with a 600 Oe field along [110]) as well 
as the minimum in specific heatP^. When increasing the 
magnetic field above 50 Oe, the magnetization starts to 
saturate, thus leading to a smaller value of M/H (See 
Fig. [2]), and the irreversibility is reduced. 

Figure [3] shows an investigation of the frequency depen- 
dence of the AC susceptibility associated with this irre- 
versibility. We measured the AC susceptibility along the 
[111] and [110] directions for samples CI and C2 respec- 
tively, and the behavior is identical for both directions. 
As seen in Figure [3j the real part x' °f the susceptibil- 
ity as a function of temperature presents a well-defined 
peak, which occurs at t he same temperature as in pre- 



vious measurements^!]. The imaginary (or dissipative) 
part x" is about 10 % of x' , indicating that a substantial 



3 




0.1 0.2 0.3 0.4 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
T(K) 

FIG. 2. (Color online) Field Cooled M/H vs. T along the 
[110] direction for sample C2 at different applied fields. Inset: 
M/H vs T with H = 5 Oe in a ZFC-FC procedure. 



amount of the sample is involved in these slow dynamics. 
For samples B et D, the peak in x' occurs in the same 
range of temperature. Nevertheless, it is broader and x" 
exhibits a more complex shape (see Suppl. Mat.). These 
features indicate that the physical origin of the freezing 
is intrinsic to the compound and that only the detailed 
dynamics might be affected by the slight differences be- 
tween the crystals. 

The AC susceptibility behavior is very different from 
a canonical spin-glass behavior. Usually, a maximum in 
x'(T) is caused by a competition between the increase 
of the magnetic polarization and the slowing down of 
the dynamics when the temperature decreases. The ex- 
istence of a distribution of relaxation times broadens the 
maximum and makes it difficult to compare %' to usual 
dynamic laws because the freezing temperature T, which 
can be defined as the x' maximuirPSl, does not provide 
a direct determination of the maximum relaxation time 
7^21. Nevertheless, in canonical spin glasses, two rela- 
tions are often used to describe the frequency depen- 
dence of the x' maximunPSl [) the Vogel-Fulcher law 
t = TQexp[E/(T ~ To)], where E is an energy barrier, 
To is a phenomenological parameter, To is the intrinsic 
relaxation time, and t is related to the measurement fre- 
quency t = l/(27r/), and ii) the dynamic scaling law 
t = 7o[(T — T c )/T c ] zl/ , where T c is the transition temper- 
ature and zv is the dynamic critical exponent. In the case 
of Tb2Ti2C>7, neither of these laws fits the temperature 
dependence of the x' maximum, thus indicating the ab- 
sence of a canonical spin-glass transition. However, the 
law r = T exp[(E/T) a ] fits the data with r 1.1 x 10~ 9 
s, E « 0.91 K and a ps 2 (See the inset of Fig. [3). Such a 
law has been proposed i n the context of a zero tempera- 
ture spin-glass transitio n 1 32 ! 33 !. In this case, the maximum 
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FIG. 3. (Color online) AC susceptibility measurements along 
the [110] direction for sample C2. Top: \' vs - T together 
with M/H vs. T in a ZFC-FC procedure with H DC = 5 Oe. 
Bottom: x" vs - T. Xac is plotted for frequencies between 
5.7 mHz and 570 Hz with H A C = 0.5 Oe. Inset: r = 1/(2tt/) 
vs. 1/T, where T is the temperature of the maximum of x' 
(red circles) and of x" (blue squares). The line is a fit to the 
equation r = r exp[(£'/T) <T ] with To = 1.1 X 10~ 9 s, E = 0.91 
K and a = 2. 



of the susceptibility is related to a glassy behavior rather 
than to a spin-glass phase transition. We can note that 
although this detailed analysis of the frequency depen- 
dence of the x' peak shows that Tb2Ti2C>7 does not ex- 
hibit a canonical spin-glass transition, the peak shift per 
decade frequency AT/[TA(log /)] is found to be 0.06- 
0.08 (depending on the frequency range) and is thus sim- 
ilar to that for insulating spin-glasses^. 

In addition to this analysis of x', it is interesting to fo- 
cus on the dissipative part x" (See the bottom of Fig.|3|. 
In a simple thermally activated process with a single en- 
ergy barrier E, the frequency dependence of the maxi- 
mum of x" as a function of 1/T would follow an Arrhc- 
nius law t = Toexp(E/T). This is not the case (See the 
inset of Fig. |3| , which suggests a more complicated pro- 
cess. The unusual shape of the x"(T) curves confirms 
it. At low frequency, / < 21 Hz, the x" onset occurs at 
low temperature (T < 0.4 K) and the peak is quite sharp 
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and asymmetric. At higher frequency, we enter a second 
regime in which the peak is much more rounded. This 
might be due to an additional relaxation process at these 
frequencies. Indeed, an intriguing feature arises in this 
regime: The value of x" decreases very slowly when in- 
creasing the temperature and falls to zero at much higher 
temperature. For example, for frequencies larger than 
570 Hz, the x" onset occurs at temperatures larger than 
4 K. This "high frequency" regime is observed in all the 
samples and indicates that slow dynamics emerges at a 
much higher temperature than the effective freezing (200 
mK) of the compound. Since it was shown that corre- 
lations at 2.5 K extend only over a single tetrahedrorP, 
this high temperature dissipation might be related to the 
dynamics of these single tetrahedra. 

These two regimes can be related to previous micro- 
scopic measurements, although in different time scales. 
When decreasing the temperature from 10 to 1 K, the 
/iSR asymmetry shows a change in shape and the /iSR 
relaxation rate 1/Xi increases rapidlySES. Meanwhile, 
the neu tron scattering quasi-elastic energy linewidth 
decrease d 23 * 27 ^ . Together with the onset of x" around 4 K 
at "high" frequencies, these results indicate that, in this 
temperature range, several characteristic times exist in 
the system (from 10~ n to 10 -2 s) and that all of them are 
increasing. At lower temperatures, below 1 K, the /iSR 
relaxation rate starts to saturate. In addition, neutron 
spin echo (NSE) measurements show that, below 400 mK, 
more than 10 % of spins are nearly frozen within the NSE 
time scale window (< 1 nsj^. The low temperature peak 
in AC susceptibility measurements might be induced by 
these "frozen" spins at the NSE scale. These observations 
show that the original magnetic ground state observed 
in Tb2Ti 2 07 results in a rich dynamic behavior over a 
very broad time range, from less than the nanosecond 
to the quasistatic regime. We can wonder whether these 
dynamics could be induced by the existence of multi- 
spin excitations. Indeed, loops with a broad distribution 
of lengths were identified in the Coulomb phase of spin 
ice^. In kagome systems, similar loops were shown to 
have temperature dependent relaxation times and to in- 
duce a dynamical freezing^. 

It is worth noting that the shape of our AC curves, 
especially the x" part, is very similar to recent mea- 
surements in Tb2Sn! 8 Ti .2O7E3. It shows that the slow 
dynamics probed by AC susceptibility is robust against 
deformation and confirms that is is not simply due to 
defects or impurities in the samples. 

Finally, the AC susceptibility measurements with a DC 
applied field show a very complex behavior (See Fig. [4} . 
i) The x' signal is much smaller than the corresponding 
M/H amplitude, showing that a part of the spins are eas- 
ily polarized by the DC field, become blocked and cease 
to respond to the AC field. In parallel, the application 
of a DC field strongly suppresses x" a t l° w frequency. 



These results are related to saturation effects and to the 
suppression of the irreversibility in the ZFC-FC curves 
reported above when increasing the magnetic field, ii) 
On the contrary, the high temperature part of x" at high 
frequency (/ > 21 Hz) remains important, even in pres- 
ence of a 1000 Oe DC field (not shown). The origin of 
this high frequency broad signal remains unclear, but it 
might be a key feature in understanding the Tb2Ti207 
exotic magnetic state, especially since it seems robust 
against magnetic field. 
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FIG. 4. (Color online) x' (left axis) and x" (right axis) vs. 
T along the [111] direction for sample CI, with Hac = 4 Oe 
and / = 11.1 Hz, in different applied DC fields. 

In conclusion, we have shown that the magnetization 
curves of Tb2Ti 2 07 measured along the [111] direction do 
not show evidence of a plateau down to 57 mK, contrary 
to the predictions of the quantum spin ice modeP^, and 
in agreement with Legl et a/P^l. Moreover, the [110] mag- 
netization curves are very similar to the [111] data. Be- 
low 250 mK, the ZFC-FC curves present an irreversibility 
which is reduced when increasing the applied magnetic 
field. Our analysis of the frequency dependence of the AC 
susceptibility leads us to conclude that these dynamics do 
not correspond to a canonical spin-glass transition, but 
rather to a freezing. The unusual features observed in the 
X" part of the susceptibility up to 4 K show that the slow- 
ing down of the dynamics occur at much higher temper- 
ature than the freezing temperature. The compatibility 
between our measurements on four samples and previous 
studies make us believe that these features are intrinsic, 
and that the coexistence of fast and slow dynamics is a 
clue in understanding the physics of Tb 2 Ti 2 07. 

We would like to thank R. Ballou, B. Canals, J. -J. 
Prejean and V. Simonet for useful discussions. 
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Supplemental Material 

In studies of Tb 2 Ti 2 0y, various specific heat responses have been reported in different sampleJ^l. We have 
measured the magnetization and AC susceptibility of several samples with different specific heat characteristics. 
This Supplemental Material presents these measurements performed on four samples, called B, CI, C2 and D. 
No important differences were observed between the samples, thus confirming our conclusions. 



I. SAMPLE PREPARATION 

Polycrystalline Tb 2 Ti 2 0y was first synthesized from 
Tb 2 3 and Ti0 2 powders (4N and 4N5). Sample C was 
made using commercial Tb 2 03. For samples B and D, 
Tb 2 03 was obtained by reduction of a commercial pow- 
der of Tb 4 7 (4N) under a flow of hydrogen and argon 
gas in the molar ratio 1/9. 

For the three crystals an initial heat treatment of the 
Tb 2 3 and Ti0 2 powders to 1200°C was followed by 
a second treatment up to 1350° C with an intermediate 
grinding and compaction so as to obtain a dense rod. 
The crystals were subsequently prepared from the rods 
of Tb 2 Ti 2 07 using the traveling floating zone technique. 
Crystals B and D were grown under one atmosphere of 
argon with growth velocities of 8 and 3 mm/h respec- 
tively, and crystal C was prepared under oxygen gas with 
a velocity of 7 mm/h. 
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FIG. 1. (Color online) M/H vs T in ZFC-FC procedure for 
sample D, measured in an applied field of 100 Oe along [111]. 

These crystals were characterized by heat capacity 
measurements^! We recall that sample D presents an 
anomaly in the specific heat around 400 mK, which is 
not associated with any anomaly in the magnetization 
at this temperature (See Figure [I]). The specific heat of 
samples B and C are very similar, but sample B has a 
slightly higher residual entropy. 



II. MAGNETIZATION MEASUREMENTS 
ALONG [111] IN DIFFERENT SAMPLES 
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FIG. 2. (Color online) First magnetization curves M vs. jj,oH 
at 80 mK along the [111] direction. Red dashed line: sample 
B. Blue dotted line: sample CI. Green line: sample D. 

Figure [2] shows the magnetization curves at very low 
temperature when the field is applied along the [111] di- 
rection for three samples: B, CI and D. We can see that 
samples B and CI present almost identical magnetiza- 
tion curves, whereas sample D shows a slightly different 
curvature around 0.3 T. To clarify these differences, we 
have plotted the derivative dM /dH as a function of field 
in Figure [3} 

III. MAGNETIZATION MEASUREMENTS 
ALONG [110] VS. [Ill] 

Figure [4] presents the magnetization curves of samples 
CI and C2 measured along the [111] and [110] directions 
respectively (Crystal C is the the crystal with the small- 
est residual entropy). 

The initial slope is very similar in both directions. At 
low temperatures, the magnetizations cross three times. 
Furthermore, the magnetization along the [110] direc- 
tion is slightly larger than along [111] at 8 T, even at 
4.2 K. The latter observation is consistent with other 
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measured magnetization curves^ but was not repro- 
duced by calculation^. 



IV. AC SUSCEPTIBILITY IN SAMPLE B 

The Zero Field Cooled - Field Cooled (ZFC-FC) curve 
of the DC susceptibility measured in sample B in a small 
field (i?DC = 5 Oe) presents an irreversibility below 250 
mK, in a very similar way as samples CI along [111] and 
C2 along [110] (See Figure [5^a)). Associated with this 
irreversibility, the real part x! °f the AC susceptibility 
presents a broad peak which depends weakly on the fre- 
quency, as shown in Figure [5ja). The imaginary part x" 
shows a more complex behavior (See Figure [5](b)) . The 
fact that the xac peaks are better defined in samples CI 
and C2 may indicate that the distribution of relaxation 
times is narrower in these samples than in sample B. 
The presence of local inhomogeneities, by modifying the 
local energy environment of the spins might extend such 
a distribution, resulting in broader features in AC sus- 
ceptibility measurements. Interestingly, for frequencies 
larger than 21 Hz, x" presents the same tail as observed 
in samples CI and C2 when the temperature is increased, 
falling to zero at much higher temperatures. 



FIG. 3. (Color online) dM/dH vs. fi H at 80 mK along the 
[111] direction in samples B, CI and D. 
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FIG. 4. (Color online) First magnetization curves M vs. fioH 
at 80 mK and 4.2 K. Solid lines: sample C2, measured along 
the [110] direction. Dotted lines: sample CI, measured along 
the [111] direction. 
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FIG. 5. (Color online) Susceptibility measurements as a func- 
tion of temperature along [111] in sample B. a) M/H measured 
in a ZFC-FC procedure with J/dc = 5 Oe together with the 
real part \ of the AC susceptibility, b) Imaginary part x" 
of the AC susceptibility, for frequencies between 21 mHz and 
211 Hz, with Hac = 0.4 Oe. 
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